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It was recently suggested that the sign of particle drift in inhomogeneous temperature or 
turbulence depends on the particle inertia: weakly inertial particles localize near minima 
of temperature or turbulence intensity (effects known as thermophoresis and turbophore¬ 
sis), while strongly inertial particles fly away from minima in an unbounded space. The 
problem of a particle near minima of turbulence intensity is related to that of two par¬ 
ticles in a random flow, so that the localization-delocalization transition in the former 
corresponds to the path-coalescence transition in the latter. The transition is signaled by 
the sign change of the Lyapunov exponent that characterizes the mean rate of particle 
approach to the minimum (which could be a wall or another particle). Here we solve 
analytically this problem for inelastic collisions and derive the phase diagram for the 
transition in the inertia-inelasticity plane. An important feature of the phase diagram is 
the region of inelastic collapse: if the restitution coefficient of particle velocity is smaller 
than some critical value, then the particle is localized for any inertia. We present direct 
numerical simulations which support the theory and in addition reveal the dependence 
of the transition of the flow correlation time, characterized by the Stokes number. 


1. Introduction 

Transport of inertial particles in an inhomogeneous turbulence or in a temperature 
gradient is of great importance for various industrial and natural processes. There is 
the widely known thermophores is: the tendency of the particles to migrate in the direc¬ 


tion of decreasing temperature (jVan Kampen Nl988t IVan Milligen et al. 1120051) . A simi¬ 


lar phenomenon for iner t ial particles in non- uniform turbulence is called turbophoresis 
(|CaDoraloni et a^.1ll975l : iReeks I Il983l l2014l) . The standard turbophoresis was consis¬ 
tently described within the gradient tra nsport model for the particle concentration un¬ 
der the assumption of local equilibrium ( Reeks 19831) . In that model, the turbophoretic 
velocity arises as a specific term in the expression for particle current which is directly 
proportional to gradient of the local turbulence intensity. The turbophoretic velocity is 
directed towards decreasing turbulence so that given zero concentration gradients and 
no net fluid flow in any direction, particles will migrate from high to low turbulence in¬ 
tensities. If turbophoresis is strong enough to overcome the effect of turbulent diffusion, 
the particles turn out to be localized near the minimum of turbulence. It was recently 
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suggested (jBelan et aZ.ll2014l ) that particle migration in an unbounded space could be ac¬ 
tually opposite- away from minima, if the particles are inertial enough. More specifically, 
the phenomenon of reverse turbophoresis takes place for the inertial particles placed in 
the vicinity of turbulence minimum provided that the particle mean free path is larger 
than distance from this minimum. In this case the local equilibrium assumption fails and 
there is no reduced description of particle transport in terms of spatial concentration. 
Thus, in contrast to the standard turbophoresis, the reverse turbophoresis is a nonlocal 
phenomenon, not associated with any local turbophoretic current. 

The change of direction of particle migration leads to localization-delocalization tran¬ 
sition upon the change of inertia. It also means separation: when time goes to infinity, 
particles with low inertia go to a minimum of turbulence and concentrate there, while 
sufficiently inertial particles escape to infinity. In many cases of interest, there is a wall, 
which corresponds to the minimum either of turbulence (in wall-bounded flows) or of 
temperature (in furnaces, combustion chambers and kerosene lanterns). It is not known 
how the boundary conditions at the wall affect the direction of the particle drift. Here we 
derive an analytical expression for the Lyapunov exponents associated with the motion 
of inertial particles near an inelastic wall. That allows us to predict the localization- 
delocalization transition upon the change in elasticity of collisions. A central result of 
our work is a phenomenon which one might call inelastic collapse: if the restitution co¬ 
efficient of particle velocity is smaller than the critical value that we determine, then the 
particle is localized near the wall for any inertia. The theoretical predictions are in a 
good agreement with the results of numerics that we carry. 

All the results are directly translated into the description of the statistics of the 
distance between two inertial particles in a spatially smooth and tempor ally random 
one-dimensional flow (see e.g. Iwilkinson fc Mehlig I 200, 'll: Bee et al. 2008). The prob¬ 
lem of relative dispersion is of particular importance for the description of distribution 
and collisions of water droplets in clouds and for the description of planet formation. 
Here localization means tha t particles tend to approach each other and create clusters 
( Wilkinson fc Mehlig I 2003 ). Our findings mean that in one-dimensional flow particles 
always create clusters when the restitution coefficient is below a threshold. That may be 
of importance for wide classes of phenomena in industry, geophysics and astrophysics, 
from clouds to planet formation. Note, however, that it will require further work to es¬ 
tablish quantitatively how significant are the effects of inelastic inter-particle collisions 
in higher-dimensional flows. 


2. General relations 


Consider the motion of a heavy particle embedded in an incompressible fluid flow 
u(r, t) near flat impenetrable wall. We introduce a reference system with the z-axis 
perpendicular to the wall and assume that the fluid occupies the region z > 0. The 
particle is assumed to be so small that the flow around it is viscous. Then, coordinate r 
and velocity v of the particle change according to 


dr{t) 

dt 


v(0: 


dv(t) u(r(t), t) — v(t) 
dt T 


( 2 . 1 ) 


where r is the particle response time (Stokes time). The boundary condition at z = 0 
is inelastically reflecting so that at every collision the particle loses a definite part of 
its wall-normal velocity v —>■ —(dv, where v = Vz and (3 is the constant coefficient of 
restitution. 

The fluid velocity field is treated as a random function of time. Its statistics is assumed 
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to be homogeneous (stationary) in time, whereas there is no homogeneity in space: typ¬ 
ical amplitude of velocity fluctuations is assumed to depend on z-coordinate due to the 
presence of wall. The joint probability density function (PDF) of the particle’s velocity 
and coordinate along the direction of flow inhomogeneity is defined as 

p(z, V, t) = {6iz- z{t))S{v - v{t))), (2.2) 

where z{t) and v{t) are the particular solutions of ( 12 . 11 ) and the averaging is over the 
statistics of random flow. 

There are three characteristic times for an inertial particle in spatially inhomogeneous 
random flow: the velocity correlation time Tq of the fluid, the particle relaxation time r, 
and the time needed by the particle to feel the flow inhomogeneity f, specified below. We 
first consider the case when Tc <C t, f, i.e. in particular, the Stokes number St = t/tc ^ 1. 
Then the fluid velocity field can be treated as short-correlated and PDF p.2p is the 
solution of the Fokker-Planck 


dtp = -vdzP + ldy{vp) + ^'^K{z)dlp , (2.3) 

where 7 = 1 /r and the effective diffusivity k{z) = {uz{z,t)uz{z,Q))dt describes the 
non-uniform intensity of turbulence. The boundary condition at the inelastic wall is 
dictated by particle number (or probability) conservation: the outcoming flux of particles 
at the boundary is balanced by incoming flux. Thus 

p{z = 0,v,t) = p{z = 0,—v/f],t) for > 0. (2.4) 

Our paper is devoted to the case of the quadratic profile of the diffusivity 

k{z) = pz"^, (2.5) 


where p measures the intensity of fluid velocity fluctuations. This model serves as a 
generic profile of the eddy diffusivity near a minima of intensity of random flow. Note, 
however, that for large-Reynolds wall-bounded turbulence, such model can have only 
qualitative use, since the effective diffu sivity behaves as z'^ in th e viscous sub-layer and 


2007 ). 


as z in the logarithmic boundary layer (iMonin fc Yaglom 

The equations (ESI) and ( 1 ^ also constitute a standard one-dimensional model to 


describe the PDF for relative motio n of two particles at viscous scale of turbulence 
( Falkovich et aZ.lE002l : iBec et q,/.II 2008I ). In this case, z is the inter-particle distance and v 
is the relative velocity, so that dz/dt = v, rdv/dt = —v + Au. At small separation z one 
obtains Au{z, t) = zs{t) where slf) is the z-independent gradient of fluid velocity. Assum¬ 
ing s{t) to be short-correlated we find Eq. (12.31) with k(z) = pz"^ and p = {s{0)s{t))dt. 

The boundary condition (12.41) takes into account inelastic inter-particle collisions. Thus, 
the results of further consideration are applicable also to the problem of relative disper¬ 
sion of inelastically colliding particles in one-dimensional random flow with zero time 
correlation. Note that in higher dimensions the long itudinal (radial) r elative motion of 


particles is coupled with the transversal dynamics (|Piterbarg l2002t iBec et al. 12007 


Zaichik et al. 1120081) . For this reason, the influence of inelastic collisions on evolution of 


inter-particle distance in higher-dimensional flows requires a separate investigation. 

The time evolution of an arbitrary initial phase-space distribution is determined by 
Eqs. (|2.3p and (12.41) . However, an exact analytic solution of this non-stationary problem 
is known only for the triv ial case of uniform diffusivity and ideally reflecting wall (/3 = 1 ). 
In ( Devenish et al. 19991) the stationary solution for constant k and arbitrary /3 was con¬ 
sidered under the assumption of a homogeneous distribution of particles in the core of 
flow. The results of their steady-state analysis indicate that inelastic boundary collisions 
lead to particle concentration near the wall. In general case of non-uniform turbulence the 
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space-dependent k determines a local time scale f(z) given by the time it takes the par¬ 
ticle flying with a rms velocity to experience a difference in average turbulence intensity. 
The time f(z) is estimated as the spatial scale I ~ k/k' of flow inhomogeneity divided 
by the mean local ’’thermal” particle velocity, which can be obtained by comparing the 
second and last terms In the rhs of (12.31) . v{z) ~ The validity of (12.31) requires 

that Tc -C f{z),T. Let us define the dimensionless measure of particle inertia as 

( 2 . 6 ) 


Equivalently, I(z) can be represented through the ratio of the particle mean free path 
and the scale of flow inhomogeneity. For k given by (j2.5l) the inertia degree is position 
independent, I = 

Our goal is to determine the long-term behaviour of inertial particles near the inelas¬ 
tic wall. Namely, we are interested in the direction of preferential particle migration. In 
the limit / <C 1 the approximation of local equilibrium is valid: the statistics of particle 
velocity is described by the Maxwell distribution with 7 k(z) playing the role of local 
temperature. Thus, the particles adjust velocity locally and do not sense the boundary 
conditions in the leading-order approximation. The particle concentration in the real 
space n{z,t) = dvp{z,v,t) on timescales t obeys the gradient transport equa¬ 
tion dtu = —dzj, where the diffusion current j is 


j = —K{z)dzn — n 


dn^z) 

dz 


(2.7) 


The second term in the right hand side is proportional to the eddy diffusivity gradient 
and represents usual turbophoretic current in the direction of decreasing turbulence. For 
quadratic ^-dependence of effective diffusivity the fluxless steady-state n{z) oc Xjz^ is 
normalizable at 0 —>■ oo so that particles concentrate near minimum of turbulence. 

For inertia degree I of order 1 and larger the particles are far from the equilibrium with 
local turbulence and gradient transpor t equation for the p article real space concentration 
is invalid. Recent theoretical analysis ( Belan et a7ll2014ll of quadratic model revealed a 
localization-delocalization transition upon the change of / for inertial particles near ide¬ 
ally reflecting wall: the direction of particle migration reverses as inertia degree is getting 
larger than some critical value ~ 0.827“^. In the next section we present a theoretical 
approach allowing us to describe this transition in terms of particle trajectories. In the 
Section |4] the approach is extended to the case of inelastic wall. 


3. Particles near the perfectly reflecting wall 

If the wall at z = 0 is elastic (perfectly reflecting), the time dependence of the particle 
position for each realization of turbulence is given by 

z(t) = z(0)exp ^p.v. y , (3.1) 


where a = vfz. The principal value (p-v.) of the integral corresponds to the moments 
of particle-wall collisions ti, for which variable a becomes singular. The representation 
dSH) Is valid for t ^ ti. 

We wish to determine the direction of particles migration. The long-time evolution of 
z is determined by the Lyapunov exponent 


A 


lim — In 

t —^-t-oo t 


z(0) 


lim —p.v. 

t —^-j-oo t 





(3.2) 
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J-i 


Figure 1. The dependence of the collision frequency on the inertia degree 1: analytical result 
(13.51) (line) and numerics for /? = 1, (o), 0.5 (■) and 0.1 (a). This plot convincingly demonstrates 
that collision rate does not depend on the coefficient of restitution ft. 


With the assumption of ergodicity, the time average converges to ensemble average, so 
that the mean value of cr in (13.211 can be calculated as (tr) = lim erP{a, t)da, where 

P{cr,t) = (S{a — (j{t))). It follows from (12.31) that probability distribution P satisfies the 
closed Fokker-Planck equation 


dtP = da [u'p] + pi^^dlp , 


(3.3) 


with U{a) = 7cr^/2 + cr^/3. Note, that dynamics of variable tr decouples from that of 2 
and V due to the quadratic z-dependence of diffusivity (12.51) . 

The cr-space has the topology of a circle with a = —oo glued to cr = +oo at z = 0. Every 
reflection from z = 0 corresponds to a jumping from = —oo to +oo. That topology allows 
the flux of probability —F going towards —oo and returning from +oo: U'P + P' = F. 
The flux independent of a corresponds to the normalizable steady state. 


P(ct) = 

i-oo 


(3.4) 


which realizes the minimum of entropy produ ction dFeig elman fc Tsveli^ 1982 ) and is 
the asymptotic state of any initial distribution ( Gawedzki et al\ 2011 ). T he constant flux 


F is the average frequency of particle-wall collisions and behaves as (jGawedzki et al. 


20111 ) 


F = 



\ /•+“ dx f 

7-00 12 



I 



(3.5) 


where Ai and Bi designate the Airy function of the first kind and second kind re¬ 
spectively. This expression comes from normalization condition for distribution (13.41) . 
= 1 . 

We can now calculate the Lyapunov exponent (j3.2l) by using steady-state PDF ()3.4I) 



aP{a)da = 




exp 



X 

4P 


X{I3 = 1J) 


(3.6) 
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This integral is negative at / < Ic~ 0.827 ^ and positive otherwise ([Wilkinson fc Mehlig 


200.‘lh . That sign change is a phase transition, called localization-delocalization transi¬ 


tion for motion of particle near a wall (when z is the distance from the wall) and path- 
coalescence transition for inter-particle dynamics (when z is an inter-particle distance). 
Note that even in the delocalized phase, despite going away exponentially, z(t) returns 
to zero with the time-independent mean rate (13.51) . 

The mean squares of the coordinate and the velocity, {z^(t)) = f z‘^p{z,v,t)dzdv and 
{v^{t)) = f v‘^p{z,v,t)dzdv, obey the identical closed equations 


d?{z'^{t)) d‘^{z‘^(t)) 2d{z^{t)) 

d^jv^jt)) ^ d‘^{v^{t)) ^ 2 d{v^{t)) 

dC ^ dt^ ^ dt 


473/3(z2(t)) = 0, 
47^/3 (t)) = 0 , 


(3.7) 

(3.8) 


which have solutions in the form + Ce^^*. The amplitudes A, B and C can 

be expressed through the initial coordinate Zq, initial velocity Vq and the eigenvalues Ai, 
A 2 and A 3 . The latter parameters are completely determined by the inertia degree I and 
provide the long-term exponential growth of both moments for any initial conditions. This 
growth may seem surprising in the case of localization; however, realizations with particles 
going away exist even in this case - those realization contribute negligibly (Inz(t)) but 
dominate (z^(t)). 


4. The case of inelastic collisions 

Assuming that at every collision with the wall the particle instantaneously loses some 
part of its wall-normal velocity we obtain v{ti + St) = —jdiviti — St) and z{ti + St) = 
Piz{ti — St), where z{ti) = 0 and St —>■ -1-0. In general, fSi should be viewed as an 
effective restitution coefficient related to energy dissipation d ue to collision-induced plas¬ 
tic deformations and particle-wall hydrodynamic interaction ( Legendre et a,l. 2006l l. The 
realistic models for these dissipation processes should involve fSi which is a function 
of the impact velocity. Note also that temporal or spatial irregularities at the surface 
of the wall may give rise to randomness of restitution coefficient. Remarkably, since 
a[ti + St) = —a{ti — St), the variable a = v/z does not feel the boundary conditions at 
the wall. This means that the probability density function P{a,t) relaxes to the same 
steady state (13.41) independently on details of particle-wall interactions. 

In this study we focus our attention on the idealized model with constant restitution 
coefficient jd. Then, it follows from the relation z{ti + St) = fdz{ti — St) that the particle 
coordinate for each particular realization of turbulence is given by 


z{t) = z(0)/3'^^*^ exp ^p.v. J a{t')dt' 


where, as before, N{t) is number of bounces up to the moment t ^ ti. 
The Lyapunov exponent can be then calculated as 


A(/3,/) = lim iln^^= lim -p.v. 


T(t')dt' 


■ In jd lim 

t—¥-\-CiO 


N{t) 


(4.1) 


(4.2) 


Since the statistics of a does not depend on jd, the first term in the rhs is just the 
Lyapunov exponent (13.61) in the problem with an ideally reflecting wall 


1 /■* 

lim -p.v. / a{t')dt' = \{\, I), 
t —>^+00 t J 


0 


(4.3) 
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P 

Figure 2. The phase diagram of the localization-delocalization transition: theory (line) and 

numerical simulations (o). 


while the second limit is the collision rate 


lim 

t—>■ + 00 


m 

t 


F{I), 


(4.4) 


which is given by (I3.5p . 

We then express the Lyapunov exponent for any inertia parameter / and the restitution 
factor 13: 

A(/3,/) = A(l,/) + F(J)ln/3. (4.5) 

Localization-delocalization phase transition corresponds to a zero Lyapunov exponent 
A(l, /) -I- F{I) In /3 = 0 which determines the transition curve (see Figure [2]) 


ln;9c(/) = 


A(/) 

F{I) 


2 




exp 




Calculating the integral for the infinite inertia parameter (/ 
zero value 


/3o = exp(--^ 


00 ) one obtains non- 


(4.6) 


We thus conclude that the curve of the phase transition in the plane P — I hits J = 00 at 
the Hnite (3 = I3q. There is no phase transition for less elastic walls where particles with 
any inertia localize. By the same token, particles coalesce if their restitution coefficient 
is small enough. This phenomenon can be called inela stic collapse in analogy to similar 
effect in driven granular matter ( Cornell et al\ 1998h . Even though beautiful formula 
(j4.6ll has been obtained within a short-correlated model, it perfectly corresponds to the 
direct numerical simulations presented in the next section for a synthetic velocity field 
with finite temporal correlations. 

Interestingly enough, the present calculations can be extended to the case of non- 
deterministic particle-wall interaction. Specifically, for a rough wall, one can treat the 
restitution coefficient of the wall-normal velocity /3 as a random variable having some 
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probability distribution function p(/3). Then, the particle coordinate evolves as 

N{t) 

n (4.7) 

i=l 

where Pi is the random restitution coefficient associated with *th bounce. We assume 
that random numbers Pi,P 2 , ■■■,Pn are statistically independent. 

From (I4T1) one obtains the following expression for the Lyapunov exponent 

, N(t) 

-^Inft. (4.8) 

2 — 1 

As before, the first term in the rhs is the Lyapunov exponent (|3.6p from the problem 
with elastic wall. The second term can be calculated as 

, N{.t) N 1 

lim -Vln^, = F(/) lim —y"\nPi = F{I) \nPp{p)dp. (4.9) 
t->+oo t ^ N^+oo N Jq 

That finally gives 

A = A(J) + F(/)(ln/3)p. (4.10) 

If (ln/3)p = ln/3p(/3)d/3 > ln/3o, the particle is localized for any inertia degree I. 

Analysis of the localization properties for the case of a velocity-dependent restitution 
coefficient remains a task for the future. 

Finally in this section, let us briefly discuss the effect of partial boundary absorption. 
We assume that on arriving at the wall the particle is absorbed with probability a and 
reflected with probability 1 — a. The wall-normal velocities of the particle just after and 
before reflection are related by u —> where restitution coefficient /3 is a random vari¬ 

able or a given function of impact velocity. We are interested in the survival probability 
(3(f) which is the probability that a particle has not yet been absorbed by wall after a 
time t. The decay rate of Q(h) is given by absorption probability a multiplied by the col¬ 
lision rate F: dQ/dt = —aF. At long timescales the probability distribution of variable 
a relaxes to steady state (13.4p providing the time-independent frequency of collisions F 
given by dSSl). This leads to the exponentially fast long-term decay Q{t) oc exp(—oFf). 
Note that the survival probability does not depend on the restitution coefficient /3 due 
to the aforementioned insensitivity of cr to the boundary condition at the wall. 


A = 


lim —p.v. 
t —>' + 00 f. 


<7{t')dt' 


lim 


z(f) = z(0) exp ( p.v. / a{t')dt' 


5. Simulations 

We conducted numerical simulations of the motion of particles subjected to one¬ 
dimensional inhomogeneous random force. In doing numerics we model the gradient 
of fluid velocity as a telegraph noise. That means that we choose u = where the 

random process u = z^{t) is a piece-wise random constant during Tc, chosen from a 
Gaussian distribution with a zero mean and unit variance. Varying we are able to 
explore finite-correlated flows as well. The equation of motion Eq. ^2TD is solved by the 
second-order Runge-Kutta algorithm with time steps much smaller than Wall reflec¬ 
tion is modeled as follows: if a particle changes the sign of z during Tc, its coordinate 
flipped to opposite while its velocity changes sign and is multiplied by the restitution 
coefficient p. To provide the relevance of the numerical simulations to our analytical 
model, we have used the the correlation time of fluid velocity Tc = 0 .002, which is much 
smaller than other time parameters of the problem, <C r, F“^,A“^. In this limit, the 
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Figure 3. The time dependence of the mean logarithm of the particle coordinate at fixed 
elasticity /3 = 1 for different degrees of inertia: (1) 1 = 1.25147 , (2) 1.21545, (3) 1.20939, (4) 
1.20332 and (5) 1.19144. 



Figure 4. The time dependence of the mean logarithm of the particle coordinate at fixed inertia 
degree I = 7.37 for different restitution coefficients: {1) 13 = 0.21 , (2) 0.205, (5) 0.20, (4) 0.195 
and (5) 0.19. 

particle displacement during Tc is small and the telegraph process reproduces the case of 
a white noise. 

We compute Lyapunov exponents to determine the direction of the particle drift 

A = lim - /in | | \ , (5-1) 

t^oo t \ z(0) / ’ ^ ^ 

where averaging is over many (typically 1000) realization of the random process. The 
Lyapunov exponent is known to be self-averaging quantity, that is having the same value 
in any given long realization as after averaging over many realizations. The localization- 
delocalization phase transition occurs when the Lyapunov exponent changes its sign. 
The hgure [T] indicates that, as expected, the collision frequency F does not depend on 
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P and is completely determined by the inertia degree I and relaxation time r. Figures [3] 
and |4] show that the mean logarithm of particle coordinate is a linear function of time. 
The slope changes its sign upon the change of inertia degree or elasticity. The phase 
diagram of localization-delocalization transition demonstrates a very good agreement of 
numerics and analytic theory, see Fig. [2j 


6. Phase transition upon the change of the Stokes number 

As stated above, there are three time parameters for an inertial particle in a random 
flow: the velocity correlation time Tc of the fluid, the particle relaxation time r, the time 
needed by the particle to feel the flow inhomogeneity f. In the previous sections the 
transition upon the change of inertia parameter I = was demonstrated in the 

limit Tc <C T, f. Here we consider Tc,t <C t and discuss the phase transition upon the 
change of the Stokes number S t = t/tc. The par ticles concentration on long timescales 
(t ^ T,Tc) obeys the equation (|Belan et a^.ll2014ll 


dtn = 


St 


1-hSt 


dl[K{z)n] -h 


1 


1-hSt 


d^[K{z)d:,n] 


( 6 . 1 ) 


For quadratic diffusivity profile k(z) = this leads to the following expression for the 
average logarithm of particle coordinate (Inz(t)) = J n(z,t)^nzdz 


dt(lnz(t)} = A = fj. 


1-St 

1-hSt 


( 6 . 2 ) 


At the critical Stokes number Stc = 1 the Lyapunov exponent vanishes and localization- 
delocalization transition occurs. The reason of this transition is that at low particle inertia 
(St < Stc) the turbophoresis can not compensate the effect of turbulent diffusion which 
tends to spread the particles throughout the fluid. On the plane x = I,y = 1/St we 
expect the phase-transition curve to start at x = Ic,y = 0 and monotonically go down 
to a; = 0, 2 / = 1. This is supported by the numerical data obtained with increasing ration 
of the flow correlation time to the Stokes time. Taking ,5 = 1, we find C = 1.209 for 
St = 500, Ic = 1.204 for St = 50 and C = 1.137 for St = 5. It is important to keep in 
mind that an accou nt of fluid incompr essibility (and multi-dimensionality) was crucial 
in deriving (16.ip by I Belan et o.l\ (I 2 OI 4 II . Therefore, the motion of particles with t < Tc 
cannot be modeled numerically with a one-dimensional fluid flow. 


7. Discussion 

We have studied analytically and numerically the relative motion of inertial particles 
in homogeneous random flow or, equivalently, the particle motion in inhomogeneous 
wall-bounded flow. The main part of our analysis was focused on the case of a short- 
correlated flow with quadratic dependence of the effective diffusivity on the distance, 
which corresponds to the inter-particle distance in the viscous interval of turbulence 
or to the motion of a single particle near a minimum without a wall. In this model, the 
direction of particle migration is determined by two dimensionless parameters: the inertia 
degree / and the restitution coefficient of particle-particle or particle-wall collisions /3. We 
have derived analytically the phase diagram for the localization-delocalization transition 
in the plane I — p. Remarkably, there is a critical value of the coefficient P below which 
particles localize for any value of inertia I. The results are confirmed by the direct 
numerical simulations of particle motion in a synthetic random flow. 

Our results do not apply directly to the near-wall region of developed turbulence. 
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where the turbulent diffusivity behaves as z* inside the viscous sub-layer and as z in 
the logarithmic boundary layer. An interesting question is whether the results provided 
by exactly solvable quadratic model can be generalized to these cases. Recent analysis 
( Belan et a/.ll20153l of inertial particles in the viscous sub-layer of wall-bounded turbu¬ 
lence (k oc z^) revealed the localization-delocalization transition at the same /3 = /3c- For 
K oc z"* we thus have the same threshold for localization at m = 2 and m = 4. Moreover, 
the critical value (14.61) is the same for inelastic collapse of trajectories under the action 
of a spati ally uniform random force, k(z) =const, as long as it is zero upon the contact, 
k( 0) — 0 ( Cornell et al. 1998ll . To address the localization properties of inertial particles 
for any m ^ 0 one needs to consider the boundary value problem (12.3112.41) , which will be 
the subject of future work. 

It is important to stress that our one-dimensional consideration with a short-correlated 
flow is a crude approximation of reality, particularly for the problem of two inertial 
particles. Further numerical and experimental work is needed to establish whether there 
is some analogue of the collapse transition leading to clustering of particles in real flows. 
This seems likely as one would expect that inelastic inter-particle collisions decrease 
the Lyapunov exponent associated with the evolution of distance between two inertial 
particles below the viscous scale of turbulence. 
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